Viscosities of the quasigluon plasma 
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Abstract. We investigate bulk and shear viscosities of the gluon plasma within relaxation time 
approximation to an effective Boltzmann-Vlasov type kinetic theory by viewing the plasma as 
describable in terms of quasigluon excitations with temperature dependent self-energies. The 
found temperature dependence of the transport coefficients agrees fairly well with available 
lattice QCD results. The impact of some details in the quasigluon dispersion relation on the 
specific shear viscosity is discussed. 



1. Introduction 

Bulk and shear rj viscosities are important properties of strongly interacting matter. Their 
firm knowledge is a necessary prerequisite for dynamical descriptions of relativistic heavy-ion 
collisions, but also in cosmology and astrophysics. The success of ideal relativistic hydrodynamic 
simulations, cf. e. g. pQ, in quantitatively describing the collective flow behaviour of the matter 
produced at the Relativistic Heavy Ion Collider led to the conclusion that at most small 
dissipative effects are present in the created medium ElSj . Later investigations by means of 
dissipative relativistic hydrodynamics, cf. e. g. El [71 IS], confirmed this observation. 

However, the quantum mechanical uncertainty principle enforces a fundamental lower bound 
on the shear viscosity [9] for any physical system. Moreover, unitarity arguments led to a lower 
limit conjecture for the specific shear viscosity [initti] ir]/s)KSS ^ h/{^7r) (Kovtun-Son-Starinets 
bound, s is the entropy density). The bulk viscosity, in contrast, is exactly zero in conformal 
theories. In QCD, however, it is only (approximately) zero at large temperatures T, while for T 
close to the deconfinement transition temperature Tc the specific bulk viscosity C/s expected 
to become large [T2 | [T3l [Ti] . 

The viscosity coefficients were recently evaluated by means of non-perturbative first-principle 
numerical lattice QCD calculations [151 [HI H?! |T8] . Other analytic approaches are based on the 
Kubo- formalism [19] or on linearized Boltzmann kinetic theory \19\ \20 \ [2T| I22j. In fact, it was 
shown that in scalar theory \19\ [20] and in hot QED [211, (22] both approaches are equivalent. 
For small QCD running coupling, perturbative QCD (pQCD) kinetic theory results at leading 
order for shear and bulk viscosities were reported in |23[ 124] and [25], respectively. 



2. Bulk and shear viscosity coefficients 

We calculate bulk and shear viscosities of tiie gluon plasma by assuming that it can be described 
by quasigluon excitations with T-dependent self-energy n(T), cf. [26\ 127) . In local thermal 
equilibrium, this picture is comprised within a quasiparticle model (QPM), where the gluon 
dispersion relation reads = \/p'^ + n(T). n(T) = T^G^(T)/2 contains the effective coupling 

G'^{T) = IGvr^/ ^11 In [A(T — T^) /Tc]^^ with parameters A and Tg. The QPM works impressively 

well for describing equilibrium lattice QCD thermodynamics, cf. |28l [291 [30] . 

By means of an effective kinetic theory |19t I20j . the model can be extended to systems 
facing small deviations from local thermal equilibrium. In fact, as T is space-time dependent in 
general, so is the dispersion relation E{x). The corresponding energy-momentum tensor reads 

T^'^(x) = / -fj:^^^f^p^{x)p'^ {x)b{x,p) + g^^^ B{Ii{x)), where b{x,p) represents the quasigluon 
distribution function entering E{x) and = {E,p). This expression for T^'^ is interwoven 
with a Boltzmann-Vlasov type kinetic equation [19j. The mean field term B{Il) is necessary 
for satisfying locally energy and momentum conservation, df^T^'^ = 0, established when [26l [27] 

dB/dH = — ^ / (27r)aE(x) ^^^'P)" This also guarantees thermodynamic self-consistency of the 
approach in thermal equilibrium and E = 6T^^/6b as a basic principle of statistical mechanics. 

Expanding h{x,p) = b^{x,p) + 6b{x,p) for small disturbances 6b from equilibrium b^, T'^^ is 
decomposed in an equilibrium part and first-order corrections thereof. Using the relaxation time 
approximation for determining Sb in the latter, one reads off [27] for r] and C in the local fluid 
rest frame by comparison with the phenomenological definitions [31] 



^ = T^/|f3-(^)[l + ^-VT)]^.-, (1) 
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where n{T) = (i(e^°/"^ — with degeneracy factor d = 16, dP/de is the squared speed of sound 
and r the relaxation time. These results generalize previous work [321 133] . in which excitations 
with constant mass M were considered. While the expression for r] looks formally the same, 
cf. also |34j . C, in Eq. ([2]) contains a combination of E^ and self-energy derivative in contrast 
to |32| , I33| , [53]. It is this particular combination which guarantees — )• in the conformal limit. 



3. Numerical results 

For numerically evaluating Eqs. ([1]) and ([2]), we employ as ansatz for the relaxation time 
T = [aiTG^(T) In(a2/G^(T))]^^, which is inspired by considerations in [35]. It turns out that 
with this ansatz for r, both r] and C resemble at large T the behaviour with temperature and 
coupling known from pQCD [23[ [2^ I25j . The parameters in G'^(T) are adjusted such that 
thermal equilibrium lattice QCD results for s/T^ from [361 137j are described by the QPM. Local 
thermal equilibrium quantities, such as s, are obtained by inserting n(T) for b^{x,p) into T^^'^{x), 
cf. [27]. Furthermore, we fit the parameters in r as oi = 2.587 x 10~^ and 02 = (//*/T)^ (Fit 1) 
or as ai = 3.85 x lO'^ and 02 = 2{^^/Tf (Fit 2), where n^/T = 2.765 as in pQCD [23]. QPM 
results for tj/s and the corresponding relaxation times are exhibited in Fig. [TJ while results for 
C/s are depicted in Fig. [21 

Good agreement with available lattice QCD results is found. In particular, r]/s exhibits a 
minimum in the vicinity of Tc, which is common for a variety of liquids and gases [111 I40| . We 
find = 0.096 at T = 1.22 Tc for Fit 2. The exact location of the minimum is driven by 

02 in T. With increasing T, r]/s increases rather mildly, while for T — ?> it dramatically rises 
in line with the relaxation time and the phenomenon of critical slowing down. 




Figure 1. Left: Specific shear viscosity as a function of T/Tc- Dash-dotted and sohd curves 
exhibit QPM results for Fit 1 and Fit 2, respectively, compared with lattice QCD data from |17j 
(full squares), [TS] (diamonds and triangles) and [15] (open squares, open and full circles). The 
dotted line exhibits the KSS bound. The grey band at large T depicts the pQCD result r/jvLL/s, 
where t/tvll is from [24] and s is from [3S]. The dashed curve (almost on top of the solid 
curve) shows rj/s for Fit 2 for constant mass M = 200 MeV instead of 11 (T) in . Right: 
Corresponding relaxation times r. The behaviour is crucially determined by the parameter a2- 

For comparison, rj/s for Fit 2 based on Eq. ([1]) but for particles with constant T- independent 
mass M = 200 MeV is also exhibited in Fig. [1] (dashed curve). Only the region around the 
minimum in rj/s is affected (the minimal value is about 14% larger than {r] / s)'^in) ^ while at 
large T details in become irrelevant. Increasing M generally leads to a decrease in rj/s. For 
M = 1 GeV a minimal value 9% smaller than ir]/s)min found. 
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Figure 2. Specific bulk viscosity 
as a function of T/Tc- Dash-dotted 
and solid curves exhibit QPM results 
for Fit 1 and Fit 2, respectively, 
compared with available lattice QCD 
data from [18] (squares) and [16] 
(circles). The dotted curve shows for 
comparison results from holographic 
QCD [39J. The grey band at large T 
depicts the pQCD result Clo/s, where 
Clo is from [25j and s is from [38j. 



4. Conclusion 

We derived expressions for bulk and shear viscosities within relaxation time approximation to an 
effective kinetic theory for the gluon plasma assumed to be composed of quasigluon excitations 
with T-dependent dispersion relation E^. Medium effects show up in t] only implicitly via 
while they are prominently evident in Specifying the relaxation time and adjusting model 
parameters to thermal equilibrium lattice QCD results, we find good agreement with available 
lattice QCD data for rj/s and C/s for T close to Tc. The infiuence of moderate constant masses 
M on rj/s is observed to be small and focussed on the region around its minimum. At large T 
our results resemble the parametric dependence on T and coupling known from pQCD. 
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